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(57) A chromosome using each of a plurality of pa- 
rameters of a physical model of a semiconductor element 
as a gene is defined and the parameters are optimized 
using a genetic algorithm based on the characteristics 
measurement data of the semiconductor element fabri- 
cated by way of trial. In the selection processing of the 
genetic algorithm, a sum of a first evaluation value based 
on linear scale data and a second evaluation value based 
on logarithmic scale data is employed as the evaluation 
value of the chromosome. 
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Description 

Technical field 

5 [0001] The present invention relates to a parameter tuning device, and in particular, relates to a parameter tuning 
device capable of tuning a large number of parameters in a physical model of a semiconductor device such as a transistor 
in a short time. 

[0002] When manufacturing an LSI. first, several samples of transistors (MOSFET) having different sizes such as a 
channel length L and channel width W of gates of the transistor (MOSFET) are manufactured on a processing line. Next, 

10 from the results of nieasurement of electrical characteristics of the test devices, a plurality of parameters in a physical 
model of the transistors is tuned so as to match to characteristics of transistors manufactured on the processing line 
with high precision. A!so, using the physical model of the transistors, simulations of various kinds of LSIs (transistors) 
manufactured on the processing line are' performed using a well-known circuit simulator such as SPICE. 
[0003] A physical model of transistors represents relationships such as Vg (gate voltage), Vd (drain voltage), and Id 

15 (drain current) with equations including variables such as a gate channel length L, a channel width W, and a plurality of 
parameters, and a large number of-models have been proposed. In the above simulations, for example, a typical and 
well-known BSIM (Berkeley Short Channel IGFET Model) is used. 

The BSIM is formed of a large number of equations, and the number of parameters to be tuned is more than 50. Regarding 
the physical models of transistors and conventional parameter tuning methods, because they are described for example 
20 in the following document, their detailed explanations are omitted.' 

Non-patent document 1 : Toru Toyabe, ed., "MOSFET Modeling and BS1M3 Users Guide", February 28, 2002. published 
by Maruzen 

[0004] Conventionally, parameter tuning devices which automatically perform parameter fitting (tuning) of physical 
models including plural parameters using genetic algorithms from experimental results, and the like, have been proposed. 
25 For example, in Patent Document 1 filed by the present inventors, a general parameter tuning device which automatically 
performs parameter tuning process of physical models including plural parameters using genetic algorithms has been 
proposed. 

Patent document 1: Japanese Patent Publication (Kokai) No. 2003-108972 
30 DISCLOSURE OF THE INVENTION 

The subject of the Invention to be Solved 

[0005] Inthe above conventional parametertuning method, becausea large numberof parameters can not be optimized 
35 simultaneously, tuning of all the parameters is perfonned by iterating an operation of first optimizing only a part of the 
parameters and then optimizing another part of the parameters with the fixed parameter. In the method, there is a problem 
of not capable of resolving optimal parameters depending on a processing order of parameters to optimize, or of requiring 
a large amount of time and labor for resolving. 

Therefore, the present inventors perfonned a study on applying the above general parametertuning process to parameter 
40 tuning of physical models of transistors. But there is a problem that tuning of parameters with efficiency and good 
precision is difficult even when applying the conventional genetic algorithm as-is to the parametertuning of physical 
models of transistors. A purpose of the present invention is to solve the above problems. 

Means for Solving the Subject 

45 

[0006] A parametertuning device of the present invention is mainly characterized in that it has means for generating 
new parameter genes by a special crossover process. Also, It has normalization means for applying to real-number 
parameters. Furthermore, it is characterized also by an evaluation means for performing evaluation of parameters so 
as to match characteristics particular to transistors (MOSFET) with high precision. 

50 

Effect of the Invention 

[0007] The parameter tuning device of the present invention with characteristics described above has an effect that 
it becomes possible to apply genetic algorithms to parameter tuning of physical models of transistors, and optimal 
55 parameters can be decided in a short time. 
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Brief Description of the Drawings 
[0008] 

Fig. 1 is a flow chart showing an overall procedure when performing simulation using a parameter tuning device of 
the present invention; 

Fig. 2 is a general flow chart showing a parameter tuning process; 
Fig. 3 is a flow chart showing a detail of a crossover process in S43; 

Fig. 4 is an explanatory drawing showing an example of a crossover method of the present invention. 

Fig. 5 is a flow chart showing a detail of an evaluation computation process; 

Figs. 6(a) and 6(b) are graphs showing Id-Vg characteristics of a transistor; 

Fig. 7 is a flow chart showing one example of an evaluation value computation in S15; 

Fig. 8 is a flow chart showing a discrete data normalization process in S17; and 

Fig. 9 is an explanatory drawing showing a method of selection of sizes of transistors to be test manufactured. 

Explanation of Reference Numerals 

[0009] 

G: Center of gravity 
P1-P3: Parent individuals 

Best Mode for Performing the Invention 

[0010] A parameter tuning device of the present invention is implemented by' creating a program for executing a 

process indicated by a flow chart to be described, and installing the program on a well-known optional computer system 
that is capable of executing it. Regarding hardware of the computer system, because it is well known, a detailed expla- 
nation is omitted. Embodiment 1 is explained below. 

Embodiment 1 

[0011] Fig. 1 is a flow chart showing an overall procecfure when perfomning simulation using a parameter tuning device 
of the present Invention. As previously noted, when manufacturing an LSI, first, in S50, several samples of transistors 
(MOSFET) having different sizes such as a channel length L and channel width W of gates are test manufactured on 
an LSI processing linie. 

[0012] Fig. 9 is an explanatory drawing showing a method of selection of sizes of transistors to be test manufactured. 
The size selection method, for example, divides between the maximum values and minimum values of L and W at equal 
inteivals, or making sides near the minimum values more fine, and selects drvisions (sizes) of the transistors to be test 
manufactured so as to become a cross shape on the L-W plane. 

[001 3] In S51 , electrical characteristics of the test-manufactured transistors are measured. Concretely, measurements 
to obtain plural sample values (discrete data) respectively concerning Id- Vd characteristic (Vb fixed). Id- Vd characteristic 
(Vg fixed), and Id-Vg characteristic (Vd fixed) are performed multiple times while changing the fixed values. 
[001 4] In S52, a parameter tuning process of a physical model of transistors is perfomned by a method to be described 
later using the parameter tuning device of the present invention so as to match the characteristics of the transistors 
manufactured on the processing line with high precision. 

[001 5] In S53, using the physical model with the parameters tuned, operating simulations of transistors having arbitrary 
channel lengths and channel widths manufactured on the processing line are performed using a well-known circuit 
simulator such as SPICE. 

By using the parameter tuning device of the present invention, high-precision parameters of the physical model can be 
obtained in a short time, and high-precision simulations can be performed. 

[0016] Fig. 2 is a general flow chart showing a parameter tuning (fitting) process using a genetic algorithm. In 841, 
one or more discrete data groups (measurement results) are prepared. In S42, N chromosomes having art^itrary pa- 
rameters of functions in a physical model of transistors as genes are prepared as an individual population. Individual 
generation is to decide the values of the genes in the chromosome and to compute an evaluation value of that chromo- 
some. N must be a number greater than or equal to 2. 

In the BSIM, although there are 50 or more parameters as previously noted, depending on the details to be simulated, 
there may be parameters fixed to default values without tuning, or may be ignored. Accordingly, the number n of pa- 
rameters to be tuned differs according to the purpose of the simulation, and although there may be 50 or more, it also 
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may be for example 1 0. Therefore, the parameters such as the number of chromosomes N or the number of generated 
children (Child) in the genetic algorithm are changed depending on the number n of the parameters to be tuned. Ac- 
cordingly, the process becomes faster as n becomes smaller. In the embodiment, for example, the number of chromo- 
somes is setto N = n X 15. In the BSIM, because a recommended range of initial parameter values is established for 

5 each parameter, an initial value is decided randomly within the recommended range of initial values for each parameter. 
[0017] In S43, chromosomes to become parent individuals are selected from the individual population generated in 
S42. The number p of parent individuals selected in the process must be less than or equal to N. Child individuals are 
generated by a crossover process to be described later from the selected parent individuals. In S44, the evaluation 
values of the child individuals generated in S43 are computed (the details are described later). 

10 [0018] In S45, the parent individuals selected in S43 and p in the order of better evaluation from the child individuals 
generated in S44 are returned to the individual population, and the rest are discarded. By this process, the chromosomes 
having low evaluation values are eliminated. In addition, a method also may be used, in which a portion of the parent 
individuals is returned as-is to the population without being subject to elimination, and a number in the amount of 
"remaining parent individuals" and the remaining parent individuals in the order of better evaluation from the child 
individuals are returned. 

[0019] In S46, it is determined whether an algorithm switching condition is satisfied. If the condition is not satisfied, it 
returns to S43, but if the condition is satisfied, it moves to S47. As the condition, whether the number of computations 
exceeds a prescribed value, or a rate of reduction of evaluation value becomes less a prescribed value can be used. In 
S47, tuning of the parameters is performed by the well-known Powell method as a local search method or otherwell-known 

20 local search method. By switching the search method in this manner, the time of parameter tuning is shortened. 

[0020] Next, a crossover process in the present invention is explained. As a conventional crossover method, a process 
of partially replacing genes of chromosomes has been performed. Although the conventional crossover method is effective 
when the genes have bit values (0 or 1 ), it is not necessarily an effective crossover method when the genes are real 
numbers. The reason.why is because a space of genotypes defined in the binary expression has a different phase 

25 structure from a space of the actual parameters after conversion to real numbers, and therefore continuity of the param- 
eters is not considered. 

[0021] For example, in the case that one bit in a binary number is inverted, the parameter afterthe inversion, depending 
on a position of the inverted bit, becomes larger when the highest'bit is inverted, but the parameter value almost does 
not change at all when the lowest bit is inverted. Because the parameters subject to tuning in the present invention are 

30 real numbers, it is difficult to perform effective tuning with the conventional crossover method. Therefore, the crossover 
method oriented toward real numbers described below is used rather than the conventional crossover method. 
[0022] Fig. 3 is a flow chart showing a detail of the crossover process in S43. The crossover method is oriented toward 
real numbers, which generates genes of child individuals from a simplex computed from genes of plural parent individuals: 
Fig. 4 is an explanatory drawing showing an example of the crossover method. In S31 , p parent individuals are randomly 

35 selected from the individual population. The value of p is desirably defined as p = n + 1 when there are n parameters to 
be tuned. In S32, the center of gravity G of the p parent individuals selected in S31 is sought. That is, an average value 
is sought for each parameter. 

[0023] In S33, a variable c is set to 1, and in S34, one child individual is generated by Equation 1 below using the 
center of gravity G and a unifomnly distributed random number. 

40 

[Equation 1] 

45 ■ ■ . 

. ^ J Qik = 0) 

r, =(u(0,l)^) 

55 

Here, p is the number of parent individuals selected, C.is a vector indicating the chromosome of the child individual 
generated, and Pk is a vector indicating the chromosome of the selected parent individual. In the present embodiment, 
the number of parent individuals selected is set to n + 1. Also. u(0,1) is a uniformly distributed random number in an 
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interval [0,1]. 



[0024] In S35, 1 is added to the variable c, and in S36. it is determined whether the variable c is greater than the 
constant "child". If the determination result is negative, it returns to S34, but if it is affirmative, the crossover process 
ends. By this processing, a "child" number of child individuals are generated. "Child" is desirably about 10 x n. 
5 [0025] Fig. 4 is an explanatory drawing showing a search range of simplex crossovers (generated range of child 
individuals) when there are 2 parameters a and P to be tuned, and the number of parent individuals randomly selected 
from the individual population is 3, The generated range of child individuals (inside of outer triangle in Fig. 4) is decided 
by multiplying by e the vectors of parent individuals up to P1-P3 from the center of gravity G, and a child individual is 



generated using a unifonn random number from that range. The recommended value of e Is ^(^p + 1) when the 

number of parent individuals is p. When the number of parameters is 3 or more, the generated range of child individuals 
becomes a space inside a hyperpolyhedron surrounded by plural hyperplanes. 

[0026] By using such a crossover method, in response to problems in which the parameters to be subject to tuning 
are real numbers, the parameters can be handled positively, and effective tuning can be performed. Being able to handle 
positively means that an individual in the vicinity of the parameter space is also in the vicinity of the gene space. Such 
a crossover method has characteristics that it is robust to the dependencies among variables and it does not depend 
on the method of scaling. It is suitable for tuning of parameters of functions in an electrical characteristics model of 
transistors in which the dependencies among the parameters are strong and there are a large number of parameters 
having different scaling. 

[0027] Usually, in the genetic algorithms, a process called mutation is performed in addition to the crossover Mutation, 
in the case of conventional genetic algorithms that handle discrete binary data, perfomns an operation of inverting a part 
of genes of a chromosome. Also, in the case of genetic algorithms that handle real numbers, an operation of adding a 
nonnal random number generated according to the normal distribution N(0, o^) to each gene of the chromosome has 
been proposed. However, because the crossover method of the present invention as noted above uses random numbers 
in the crossover process, it combines also the property of mutation. Therefore, the mutation is not performed in the case 
of using the crossover method as noted above. 

[0028] Next, the evaluation value executed in S44 is explained. The evaluation value of a chromosome is computed 
from prepared discrete data groups and estimation data groups computed as a function in an electrical characteristics 
model of transistors using genes of the chromosome as model parameters. The evaluation value is a value indicating 
how close the genes in the chromosome are to Ideal values as model parameters. 

[0029] Fig. 5 is a flow chart showing a detail of the evaluation computation process. In SI 2, the gene information of 
the chromosome is read in, and it is taken as model parameters in functions of the electrical characteristics model of 
transistors. In SI 3-1 ... Si 3-d, the discrete data groups are read in. d is the number of discrete data groups read in. 
[0030] In 314-1 ... S14-d, estimation data groups for predicting the discrete data groups are computed based on the 
model parameters input in SI 2. The estimation data groups are present in the same number as the discrete data groups. 
In SI 5, the evaluation value of the chromosome is computed from the discrete data groups and the estimation data 
groups. As the evaluation function, the square error to be described later, or the like, is used. 

[0031] Fig. 6(a) is a graph of linear scale and Fig. 6(b) is a graph of log (logarithmic) scale showing the Id-Vg charac- 
teristics of a transistor. As in the parts surrounded by the ellipses in Fig. 6, there is a part in the discrete data group, 
where it appears as almost 0 when viewed on the linear scale in Fig. 6 (a) but the values are greatly changed when 
viewed on the log scale in Fig. 6(b). Such a property is called sub-threshold characteristic. 

[0032] Because the absolute value of this part is small compared with the other parts, the absolute value of the error 
also is small, and it is difficult to optimize this part using only the normal linear scale data group. Also, although the 
sub-threshold characteristic can be optimized if the optimization is performed by preparing only a log scale data group 
in order to tune this part, the error of the other parts becomes greater and discrepancy occurs. 

[0033] Therefore, in the present invention, attention is on the fact that the electrical characteristics model parameters 
of MOB transistors are independent between the parameters for estimating the sub-threshold characteristic and the 
parameters for estimating the other parts. By the scaling process shown below, the log scale data group and the linear 
scale data group are read in simultaneously, and all the characteristics are tuned simultaneously. 
[0034] In the case when the scaling among discrete data groups is different, if the square error is taken, the data group 
having small scale comes to have smaller influence on the evaluation value. Therefore, there is a concern that the tuning 
precision may fall even when the scaling process previously described is perfomied. Therefore, in the present invention, 
the tuning precision is improved by nornializing the discrete data in each data group and unifying the scale. 
[0035] Fig. 7 is a flow chart showing one example of computation of evaluation value, in which the above scaling and 
discrete data nonnalization method are packaged in S15. In SI 6, logarithmic data converted the discrete data groups 
and the con-esponding estimation data groups to log scale are generated. !n SI 7, all the data groups are nonrialized. 
[0036] Fig. 8 is a flow chart showing one example of the normalization processing of discrete data performed in 81 7. 
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10 



In S21, the data groups are read in. In S22, conversion of the data indicated below is performed. That is, first, the 
maximum value fmax and minimum value fmin in the data groups are sought. Next, all the discrete data f(l) In the data 
groups are converted to normalized data g(i) by Equation 2 below. 

[Equation 2] 

2(1) = ^^^^'^ ^^^^ 
f 

J max J nun 



Here, g(i) is normalized data, f(i) is discrete data, fmax is a maximum value in data groups, and fmin is a minimum value 

in the data groups. By this computation, the discrete data can be normalized to within the range of [0, 1]. 

[0037] Returning to Fig. 7, in S1 8, the evaluation value of only the linear scale data groups is computed and is taken 

15 as A. A is a sum of the square errors between the normalized discrete data and the normalized estimation data. In S1 9. 
the evaluation value of only the log scale data groups is computed and is taken as B. B is a sum of the square errors 
between the normalized logarithmic discrete data and the normalized logarithmic estimation data. In S20, A + B is taken 
as an evaluation value of the chromosome. Although the square error is used for computation of the evaluation value, 
it also may be made such that the error ratio is sought instead of the square error. 

20 [0038] By the processing as described above, high-precision parameter tuning is possible in a short time. Also, because 
high-precision circuit simulation can be performed without test manufacturing by using the parameters in question for 
the physical model, the efficiency of manufacturing of semiconductor devices is Improved. 

[0039] Embodiment 1 is explained above, but in the parameter tuning device of the present invention, modified ex- 
amples as below also can be imagined. The BSIM is given as an example of the physical model of transistors, but the 
25 present invention can use any well-known physical model of transistors in addition to the BSIM. Furthemiore, the present 
invention can be applied to any well-known semiconductor device other than transistors. 

POSSIBILITY OF USING ABILITY IN INDUSTRY 

30 [0040] The present invention can be utilized in thesimulation of various kinds of LS Is manufactured on an LSI processing 
line using the physical models of transistors. 



Claims 

35 

1 . A parameter tuning device, comprising: 



parameter* tuning means for defining chromosomes having a plurality of parameters of a physical model of a 
semiconductor device as a gene and optimizing the parameters using a genetic algorithm based on characteristic 
measurement data of a test-manufactured semiconductor device. 



2. A parameter tuning device according to claim 1, wherein said parameter tuning means includes generation range 
deciding means for obtaining a center of gravity of a parent chromosome group in a vector space to determine a 
generation range of a child chromosome group inside a hyper-polyhedron on a vector space detemiined from the 

45 center of gravity and values of the parent chromosome group in a crossover process in the genetic algorithm. 

3. A parameter tuning device according to claim 1 , wherein said parameter tuning means includes an evaluation value 
computing means for obtaining a first evaluation value based on a linear scale and a second evaluation value based 
on a log scale to make a sum of the first evaluation value and the second evaluation value as an evaluation value 

50 of the chromosome in a selection process in the genetic algorithm. 

4. A parameter tuning device'according to claim 1 , wherein said parameter tuning means includes normalization for 
unifying a scale of data in the selection process in the genetic algorithm. 

55 5. A parameter tuning device according to claim 1 , wherein said parameter tuning means includes search method 
switching means for switching to local searching means when a parameter tuning process in the genetic algorithm 
satisfies a predetermined condition. 
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